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Abstract 



We present a study on the spectra and elliptic flow V2 for heavy flavor (charm 
^ ' and bottom) decayed electrons provided the relative contributions of charm and 

CN . bottom hadrons from the PYTHIA calculations. We made a simultaneous fit to both 

measured non-photonic electron spectra and V2 distributions. The results suggest 
that the bottom contribution is not dominant for electron pT < 5 GeV/c in the 200 
_ GeV Au+Au collisions. 
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Heavy quarks provide a unique tool to probe the partonic matter created in 
^ I relativistic heavy-ion collisions. In Au+Au collisions at RHIC, a suppression of 

I non-photonic electrons spectra [1,2] and a non-zero non-photonic electron V2 

decreasing at > 2 GeV/c [3] were observed. In the measured non-photonic 
electron spectra and V2 distributions both contributions from charm-hadrons 
and bottom-hadrons are mixed together. There are large uncertainties in the 
model predictions for charm and bottom production in high-energy nuclear 
collisions [4,5]. In order to understand the heavy quark energy loss and flow, 
which are sensitive to the properties of the hot/dense matter created at RHIC 
energy, it is necessary to separate the charm from bottom contributions. How- 
ever, the isolation of charm and bottom contributions with the semi-leptonic 
decay measurements alone remains a challenge. 

In this study, we assume the relative contributions of charm and bottom from 
PYTHIA. We then fit the measured electron spectra (Raa) and V2 simul- 
taneously. We find that the contribution of bottom is not dominant for pt 
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(electron) < 5 GeV/c. This electron px is corresponding to both D- meson and 
B-meson's pt up to 5-7 GeV/c. This simulation will provide a useful/necessary 
constraint for heavy quark production model calculations. 
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Fig. 1. Left panel: Charm/bottom-hadron and their decayed electron spectra from 
simulation fit to experiment data in p+p and Au+Au collisions. Right panel: Spectra 
ratio bands are according to da^^ /da^^ = 0.002-0.006. 

Since vector-meson decay form factor is used in PYTHIA, the spectra from 
default PYTHIA parameters are soft [6]. The heavy quark fragmentation func- 
tion is modified from the default Peterson function and the other parameter 
PARP(67) is also tuned to describe the measured non-photonic electron spec- 
tra [7]. The calculated charm-hadron spectrum is normalized to dN/dy = 
0.028 measured in p+p collisions [8] shown in the left panel of Fig. 1. A to- 
tal cross-section per nucleon-nucleon collision ratio of da^^/da^^'^ = 0.002 
is assumed to normalize the bottom-hadron spectrum. Their decayed elec- 
tron spectra are shown as dashed and dotted curves respectively. The total 
charm/bottom decayed electron spectrum shown as thick solid curve is used 
to fit the non-photonic electron spectrum from STAR EMC in p+p colli- 
sions (open squares) [1] by variating the bottom over charm cross-section 
ratio, c — >■ e from measured D° spectrum scale to total charm (dot-dashed 
curve) via form factor decay [6] and measured TOF electron spectrum (open 
circles) [9] in Au+Au 0-80% are scaled by A^bin = 293. From the fit to the 
spectra, we can see that the agreement between simulation and measurement 
is achieved when bottom contribution is small. The right panel of Fig. 1 shows 
the spectra ratio of bottom-hadron over charm-hadron (green band) and their 
decayed electron (yellow band) as well. FONLL calculation predicts the ratio 
from 0.0018 to 0.026 due to large uncertainties of the cross-section [4]. Since 
non-photonic electron spectrum is strongly suppressed in central Au+Au col- 
lions compared the spectrum in p+p collisions, scaled by number of binary 
collisions, bottom contribution seems small in Au+Au collisions. So we as- 
sume da^^/da^^'^ from 0.002 to 0.006 covered in the ratio bands. If charm 
is enhanced in RHIC energy, this ratio will be smaller [9]. With the ratio of 
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0.002, we can see at electron ~ 5 GeV/c, the bottom decayed electron 
contribution is less than 20% and decrease fast when pt goes lower. 

In order to study on the bottom contribution to the non-photonic electron 
V2 in Au+Au collisions, we propose two assumptions. Firstly, we assume 
that da^^ / da^g'^ = 0.002 — 0.006. The spectra ratios of charm-hadrons over 
bottom-hadrons are the same in p+p and Au+Au collisions at RHIC. Charm- 
hadron f 2 is the same as the light hadron. Bottom-hadron V2 is the same as 
the light hadron, or bottom-hadron V2 is zero. 




Fig. 2. Panel (a): Light flavor meson V2 and baryon V2', Panel (b)(c): Assumed 
charm/bottom-hadron V2 and their decayed electron V2', Panel (d): Comparison 
between the total decayed electron V2 and the PHENIX non-photonic electron V2- 

Panel (a) of Fig. 2 shows the meson V2 and baryon V2 from Ref. [10]. The 
assumed charm-hadron V2 and bottom-hadron V2 are shown in panel (b). Their 
decayed electron V2 are shown in panel (c). The obvious mass effect of the 
bottom decayed electron f 2 pushing to high pt is seen. The total electron f 2 
from charm and bottom decay can be obtained by applying the spectra ratio 
shown in the right panel of Fig. 1. If bottom-hadron V2 is assumed as the 
dashed curve shown in panel (b), the total decayed electron V2 will be as the 
diamonds shown in the panel (d), which becomes flat at higher pt deviating 
from the PHENIX non-photonic electron V2 (green triangles) [3]. If bottom- 
hadron V2 is zero, the total decayed electron V2 will decrease as a function of 
Pt, shown as the yehow band, due to the bottom contribution according to 
the spectra ratio. The open squares are related to the da^'^ / da^^^ = 0.006. 
Therefore, with this assumption, the decreasing of non-photonic electron V2 is 
probably due to bottom contribution, but bottom does not flow. 

The second assumption is that bottom does not contribute. A smaller and 
decreasing charm-hadron V2 is assumed shown as the dashed curve in Fig. 3 at 



3 



^ 25 
20 
15 
10 



_ 1 1 1 1 


II 

charm hat 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 
Iron 




— o 


c ^ e • ■ - 

PHENIX non-photonic e : 


: 




- Baryon ~ 










Meson 














: , , , ' 




, , , , 1 , , , , f , , , , 1 , , , , 1 ; 



(GeV/c) 

Fig. 3. The decayed electron V2 from assumed decreasing charm-hadron V2 fits to 
PHENIX non-photonic electron V2 distributions. 

higher pt- Its decayed electron t>2 (open squares) is used to fit the experimental 
data (green triangles) by variating charm-hadrom V2 distribution. If charm- 
hadron V2 is smaller than light hadrons and decreasing at pt> 2 GeV/c, its 
decayed electron V2 can describe the experimental data. 

In summary, charm/bottom and their decayed electron spectra and V2 have 
been studied using PYTHIA simulation. If bottom contribution is not domi- 
nant at Pt (electron) up to 5 GeV/c, both the measured non-photonic electron 
spectra (Raa) and V2 can be described. This result should provide some con- 
strain to model calculations on heavy flavor productions though we reached 
this conclusion with two hypotheses. The final solution to the problem has to 
come from the direct meansurement of reconstructed charm-hadron distribu- 
tions. 
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